Corrections to the universal behavior of the Coulomb-blockade peak splitting for 
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Building upon earlier work on the relation between the dimensionless interdot channel conduc- 
tance g and the fractional Coulomb-blockade peak splitting / for two electrostatically equivalent 
dots, we calculate the leading correction that results from an interdot tunneling barrier that is not a 
delta- function but, rather, has a finite height Vb and a nonzero width £ and can be approximated as 
parabolic near its peak. We develop a new treatment of the problem for g <C 1 that starts from the 
single-particle eigenstates for the full coupled-dot system. The finiteness of the barrier leads to a 
small upward shift of the /-versus- g curve for j< 1. The shift is a consequence of the fact that the 
tunneling matrix elements vary exponentially with the energies of the states connected. Therefore, 
when g is small, it can pay to tunnel to intermediate states with single-particle energies above the 
barrier height Vo- For a parabolic barrier, the energy scale for the variation in the tunneling matrix 
elements is huj, where uj, which is proportional to \/^o/C> I s the harmonic oscillator frequency of 
the inverted parabolic well. The size of the correction to previous zero-width (£ = 0) calculations 
depends strongly on the ratio between fkv and the energy cost U associated with moving electrons 
between the dots. In the limit g — *■ 0, the finite-width /-versus-*; curve behaves like (U/Hu>)/\ \ng\. 
The correction to the zero-width behavior does not affect agreement with recent experiments in 
which 2nU/hui ~ 1 but may be important in future experiments. 

PACS: 73.23.Hk,73.20.Dx,71.45.-d,73.40.Gk 
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I. INTRODUCTION 

The opening of tunneling channels between two quan- 
tum dots leads to a transition from a Coulomb block- 
ade characteristic of isolated dots to one characteristic of 
a single large composite dot.El For a pair of electrostat- 
ically identical quantum dots, this transformation can 
be chronicled by tracking the splitting of the Coulomb 
blockade conductance peaks as a function of the conduc- 
tance through the interdot tunneling channels. If one 
assumes a single common value for the conductance in 
each tunneling channel (an assumption that is exactly 
fulfilled for a spin-symmetric system of only two chan- 
nels, one for spin-up electrons and the other for spin- 
down electrons), one can divide the peak splitting by its 
saturation value and loak_for the relation between two 
dimensionless quantities:Ln3 the fractional peak splitting 
/ and the dimensionless channel conductance g .13 

For g <C 1, previous theoretical workOBO'B has treated 
the coupled-dot problem via a "transfer-Hamiltonian ap- 
proach, "E!J in which states localized on one dot are con- 
nected to those localized on the other by hopping matrix 
elements. (Here localized signifies that a state is entirely 
restricted to one of the two dots.) The hopping matrix 
elements have been treated as constant, independent of 
the states connected, as they would be if the interdot bar- 
rier were a delta-function potential, having infinite height 
and zero width. For such a barrier, the leading small-g 
behavior of the fractional peak splitting is given by 
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(1) 



where N c ^ is the number of separate tunneling modes 
(spin-up and spin-down channels are counted separately) . 

The superscript of /i^_ tells us that this is the leading 
term in the weak-coupling limit. The subscript further 
specifies that this term is calculated for a tunneling bar- 
rier of effectively zero width (£ = 0) and therefore, by 
implication, of infinite height. For g > 0.2, subleading 
terms, which are higher-order in g, contribute signifi- 
cantly to the zero-width peak splitting, and the first set 
of these, which consists of terms proportional to g 2 , has 
been calculated in previous work.tj 

In this paper, we calculate a different correction to the 
£ = 0, first-order in g result which arises from the fact 
that a realistic barrier possesses a finite height Vb and a 
nonzero width £. For such a realistic barrier, the hop- 
ping matrix elements that move electrons between the 
dots are not independent of the states they connect, and, 
for small g, they depend exponentially on the energies of 
the states. As a result of this exponential dependence, 
in the weak-coupling (§ < 1) limit, it can pay to tun- 
nel to intermediate states with energies above the bar- 
rier, and the leading term in the fractional peak splitting 
then behaves as (U/W)/\\ng\, where U is the interdot 
charging energy, which measures the capacitive energy 
cost of moving electrons between the dots, and W is the 
characteristic energy scale over which the hopping matrix 
elements change from their values at the Fermi energy 

We examine specifically the case of a finite-width in- 
terdot barrier that can be treated as parabolic near its 
peak. We find that, for such a barrier, the energy scale 
W is equal to hwjl-K, where u> is the harmonic oscillator 
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FIG. 1. (a) Schematic diagram for a single orbital-mode 
connection between the two dots. Over a distance of order 
£, the connection narrows to a minimum width on the or- 
der of the Fermi wavelength Af. (b) "Box- like" double-dot 
system with a central barrier. Hard confining walls are lo- 
cated at a distance Ldot from the barrier. The barrier is 
characterized by its height Vo and half width £. 

frequency of the inverted parabolic well. This frequency 
is proportional to the square root of the barrier curva- 
ture, which is proportional to Vb/£ 2 . It follows that the 
limit £ — ► corresponds to the limit U/W — > 0. For 
£ 7^ 0, on the other hand, it is not always true that 
U/W < 1, In fad— in recent experiments by Waugh et 
al.a, Crouch et alSi, and Livermore et al.a, it appears 
that U/W is roughly 1. Under such circumstances, we 
find that, for a given small value of the channel conduc- 
tance (g <C 1), the fractional peak splitting / is larger 
than the zero- width splitting, /£=o, by a small but no- 
ticeable amount, and, in the extreme limit of g — > 0, the 
ratio of the finite-width peak splitting to the previously 
calculated zero-width peak splitting becomes very large. 
For intermediate values of g, on the other hand, the pri- 
mary effect is a small increase in / accompanied by a 
reduction of the slope of the /-versus-g curve. 

To find the leading term in the finite-width fractional 
peak splitting we adopt a stationary-state approach!^ in 
which the first step is to solve for the single-particle eigen- 
states of non-interacting electrons moving in the electro- 
static potential of the coupled dots. The capacitive inter- 
actions between the electrons are then expressed in terms 
of these non-interacting double-dot eigenstates, and the 
off-diagonal elements of the interactions are treated per- 
turbatively. The leading term in the finite-width frac- 
tional peak splitting, is determined by finding the 
value for p = 1 of a more general quantity f^ x \p), where 
p is a dimensionless parameter (defined by Eq. || below) 
which is a measure of the bias asymmetry between the 



dotsfii In the limit U/W — > 0, the zero-width result, 
/i=0> i s recovered. For finite U/W, an approximate an- 
alytic calculation demonstrates the limiting l/|ln<7| be- 
havior, which is confirmed by numerical results. For the 
particular choice U/W — 1, which corresponds to recent 
experimentsJjQ as well as for various other choices of 
the ratio U/W, the leading term in the fractional peak 
splitting is computed numerically as a function of g. It 
is confirmed that the condition £ ^ leads to an up- 
ward shift of the peak splitting for weakly coupled dots 
(g As g becomes larger, the effect of allowing £ ^ 

becomes less dramatic, and, for U/W ~ 1, previous pre- 
dictions for the fractional peak splitting at intermediate 
values of g are essentially unaltered. 

The structure of this paper is as follows. Sec. II de- 
velops the stationary-state approach for calculating f(p). 
Sec. Ill implements this approach for a parabolic interdot 
barrier, verifying the 1/| lng| behavior of the g — > peak 
splitting that arises for £ =/= and putting the finite- width 
calculation in the context of earlier work. Sec. IV sum- 
marizes the results and comments on the possible effects 
of £ ^ when the dots are strongly coupled (g ~ I). 



II. THE STATIONARY-STATE APPROACH 

In order to solve for f(p) via the stationary-state ap- 
proach, we make the problem one-dimensional by con- 
sidering a smooth, adiabatic interdot connection [see 
Fig. 1(a)] which, for simplicity, we presume to contain 
only one transverse orbitaL-mode that lies near or be- 
low the Fermi energy i/p-llil (The use of one orbital 
mode corresponds to the spin-symmetric rN c h — 2 ex- 
periments of Waugh et aln, Crouch et alB, and Liver- 
more et alu) For such a single orbital-mode connection, 
the only parts of an electron wavefunction that can pass 
from dot to dot are those that overlap with the low- 
est transverse mode. Hence, in investigating the effect 
of the interdot connection, we can ignore all electrons 
but those in this lowest mode. We are left with a one- 
dimensional problem in which a representative electron 
with single-particle energy E moves in an effective po- 
tential V(x) = E tl -(x) + V c \(x), where E tT (x) is the spa- 
tially dependent energy of the lowest transverse mode 
and V^i(x) is the spatially dependent electrostatic energy. 
The characteristic length scale for the spatial variation of 
the effective potential is the barrier width £. 

After adding hard boundaries at a distance I/dot from 
the barrier, we have a "box-like" double-dot system with 
a smoothly varying longitudinal potential [see Fig. 1(b)]. 
The Hamiltonian consists of two components. The first, 
Hq, is a diagonal term giving the energies of the non- 
interacting, single-particle eigenstates, which form a dis- 
crete spectrum with an average level spacing proportional 
to Hv-p/Ldot near the Fermi surface, where vf is the Fermi 
velocity. The second, Hq, gives the capacitive energy 
cost of moving electrons from one side of the barrier to 
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the other 
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H c = U(n - p/2f 



(2) 



Here, h counts the electrons transferred from dot 1 to 
dot 2 (assuming, for convenience, an even total num- 
ber of electrons initially divided equally between the two 
dots). The dimensionless parameter p is a measure of the 
capacitively weighted bias and favors occupation of dot 
2 when p > O.^iNote that U equals the quantity U2 oi 
previous work.Bfl) 

We need to solve for the dimensionless channel con- 
ductance g and the fractional peak splitting /. For the 
non-interacting electrons characterized by Hq, the dimen- 
sionless channel conductance g is simply the transmission 
probability for aj-aarticle incident on the barrier at the 
Fermi energy £"p£3 / is not so easily determined because, 
in its evaluation, He is relevant. Thus, we must develop 
a means of dealing with h, which is not diagonal in the 
basis of non-interacting single-particle eigenstates that is 
the cornerstone of our approach. 

Our strategy is to switch to a basis that is simply re- 
lated to the eigenstate basis but which renders h nearly 
diagonal at energies that are low compared to the bar- 
rier. We use the fact that, for a bound system containing 
two equal potential minima, the eigenstates come in well- 
defined, discrete pairs£3 The states in these pairs have 
similar energies but opposite parities, the even-parity 
state having a lower energy than the odd-parity state. 
Thus, the non-interacting part of the Hamiltonian can 
be written in the form 



a -j 



a-3 



(3) 



where S and A are the even and odd parity indices, j is 
the pair index, and a is the spin index (which might be 
more generally regarded as a channel index). 

At lower and lower energies relative to the barrier, the 
splitting within the pairs, \Ea{]) — Es(j)\, approaches 
zero, but the spacing between pairs, \Eg(j + 1) — Ea(J)\, 
remains approximately equal to S, where 8 — Trhv-p/Ldot 
(assuming we do not stray too far from the Fermi sur- 
face). It follows that, at low energies, one can form dou- 
blets of quasi-localized states — states that lie mostly on 
one of the two sides of the central barrier — from linear 
combinations of the symmetric pad antisymmetric com- 
ponents of each eigenstate pair .113 If <j>sj{x) and (f>Aj(x) 
are the symmetric and antisymmetric eigenfunctions of 
the jth lowest-energy pair (with appropriately chosen 
overall phases), the recipe for the quasi- localized wave- 
functions is 



n ia (a;) = [<t> Sj (x) + {-iy 



(*)] 



(4) 



where the dot index a signifies that flj a (x) is primarily 
localized on the dot 1 side of the barrier if a = 1 and on 
the dot 2 side if a = 2. 



At high energies relative to the barrier, we cannot form 
quasi-localized wavefunctions from combinations of just 
two states. Nevertheless, we continue to form the linear 
combinations analogous to those of Eq. [l| We refer to 
the full set of states VLj a (x) as semi-localized to indicate 
that these states are sometimes quasi- localized (i.e., when 
they lie at low energies relative to the interdot barrier) 
and sometimes not. 

The semi-localized states constitute the complete and 
orthogonal basis that we need to render n nearly diagonal 
at low energies. Their simple relation to the double-dot 
eigenstates translates into an equally simple relation be- 
tween the corresponding creation and annihilation opera- 
tors. The semi-localized annihilation operators are given 
by 



a jaa — — J= [CSja + (-l) a+1 CAja] , 

and the corresponding expression for Hq is 



(5) 



H = ^ E{j)a) aa a 



- ^2t(j)(a] 2a a jla +H.c.) 



(6) 



where E(j) is the average energy of the pair and t(j) is 
half the energy difference within the pair: 



m = 
m = 



E A (j) + E s (J) 
2 

E A (j)-E s (j) 



(7) 



It is important to note that, whereas E(j) is in general on 
the order of the Fermi energy, t(j) is no greater than the 
average level spacing 6 = nhv-p / L^ot and becomes van- 
ishingly small in the large-dot limit (Ldot ~~ * 00 )■ The 
minuteness of t(j) will permit us to ignore it in calculat- 
ing the leading contribution to the fractional peak split- 
ting. 

We now write n in terms of the semi-localized opera- 
tors. If dot 1 corresponds to the x < side of the barrier 
and dot 2 corresponds to the x > side, we have 



dx[Q{x) - 6(-x)]V t (x)ip(x), 



(8) 



where ip(x) is the position operator and 0(x) is the Heav- 
iside step function]^ 

After writing tp{x) in terms of the semi-localized opera- 
tors aj aa , we see that h — fiQ+Snc+SfiT, where no corre- 
sponds to the h we would have if the semi-localized states 
were truly localized, Snc is the part of 8fi = (h — no) that 
does not transfer electrons from dot 1 to dot 2, and SfiT 
is the part of 5n that does effect such a transfer: 
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B{j2,ot; jx,a) 
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Sn T = ^2 5 (j2,a;ji,a)a] 2a(T a ilQC 

Here, a means "not a" and 
1 



(9) 



B(j 2 ,a 2 ;j 1 ,a 1 ) = ^J dx Vs A (*) ^ (x) 

+ (-ir +1 ct>* Aj2 (x)<t>sM]- (10) 

We have obtained the desired "semi-diagonal" form of 
n. Using 8h = dnc+Sh^ and assuming that g is small, we 
express the Hamiltonian in terms of one non-perturbative 
piece, Hq, and two perturbative pieces, and H' c : 

H o=J2 E V) a ia*<ha* + U(n - P /2)\ 

H T = - 5Z*(i)(«i2^1a + H.C.), 

H' c = U(ho - p/2)5h + U5h(h a - p/2) + U(dh) 2 . (11) 

As in Refs. 5 and 6, the fractional peak splitting is 
determined from f(p), where 



hp) 



A(Q) - Hp) 

U/A 



(12) 



and A(p) is the energy shift of the ground state of Hq due 
to the perturbations Ht^ and Hq for the given value of p, 
where < p < 1 and the total number of particles in the 

double-dot system is even. The quantity lim p ^i f(p) 

equals the fractional peak splitting /. 

Eq. [l2] tells us that we are only interested in relative 
energy shifts. Consequently, we can ignore terms such as 
(0\U(Sh) 2 \0) that are independent of p. (Here the brack- 
ets indicate an expectation value taken in the ground 
state of Hq.) Another set of irrevelant terms are those of 
the form (0\U(fio — p/2)Sn\0), which are zero due to the 
symmetry of the ground state with respect to interchange 
of the two dots. Finally, terms that contain H^ are also 
negligible because t(j) goes to zero with the reciprocal 
of the system size and, unlike Sn, H^ only connects each 
state to one other, rather than connecting each state to 
a manifold of others (see Ref. 5 for a similar situation 
with regard to odd orders in the transfer-Hamiltonian 
perturbation theory). After the above terms are omit- 
ted, it is apparent that the leading perturbative energy 
shift comes from the term that is second order in H' c . To 
lowest order in dn, this term is 



5nP 



(no - p) 2 

Kip) 



H 



P Sn 



0), (13) 



where E' (p) is the energy of the ground state of Hq and 
where Pq is the operator that projects out the unper- 
turbed ground state. A^ 2 ^(p) can easily be seen to con- 
sist of two distinct parts: a term second-order in Snc, 



which involves hopping between states semi-localized on 
the same dot, and a term second-order in Snx, which 
involves hopping between states on different dots. 

With Eq. 13, we have completed our tour of how to use 



the stationary-state approach to find both the interdot 
channel conductance g and the fractional peak splitting 
/. In order to progress further, we must adopt a model 
for the barrier that gives the en ergy dependence of the 
elements of 8n (recall Eqs. [| and 10). 



III. PEAK SPLITTING AND CONDUCTANCE 
FOR A PARABOLIC BARRIER 

We assume that the barrier in the interdot tunneling 
channel can be reasonably modeled by a parabolic one. 
For an energy barrier with peak height Voj such a model 
is plausible when Vb j^Ep » U, which is the regime of 
experimental interest .Ej The formula for a parabolic po- 
tential V(x) centered at the origin with half width £ is 
the following: 



V(x) = 




if |x| < V2£ 
otherwise. 



(14) 



A crucial energy scale for this barrier is the harmonic os- 
cillator frequency ui of the inverted parabolic well. This 
frequency is given by the formula 



huj 



1 \ {2irhf 



nV2j 2m\v£ ' 



(15) 



where 2tt/Xv = ^j2raVo/h 2 and m is the effective mass 
of the electron. 

The problem of transmission through and reflection 
from a jaaxabolic barrier is well known and exactly 
solvable .E§E3 The solutions are parabolic cylinder func- 
tions ,E3 Atid. the dimensionless channel conductance is 
given byBEfl 



1 



9 = 



1 + e -2iry(S F ) 

where E-p is the Fermi energy and 

E-V 



y(E) 



From these equations, it follows that 

(V -E F ) 1 fl-g 
= ^— In 



fko 2tt 



and, for g <^ 1, 



V 2tt 2 V2 \U ' ' 



(16) 



(17) 



(18) 



(19) 
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Eq. tells us that, even for experimental systemsi'i in 
which £ is quite small (£ ~ Af), Ep is close to Vq for 
| In <7 1 <C 27r 2 \/2- Thus, the assumption of a parabolic 
barrier appears reasonable for any measurable value of 
the interdot conductance. 

We now consider the sizes of the various energies that 
appear in our peak splitting calculations. Equation ^ 
indicates that the energy scale W for the variation of 
transmission probabilities is hw /2tt. Recalling our dis- 
cussion in Sec. I, we have 



cj> Pj (x) 



W 



2nU 



(20) 



As observed in previous work,tffl for symmetric dots, U 
equals e 2 /(Cs 4- 2Ci nt ), where Cs is the total capaci- 
tance of one of the two dots and Ci nt is the interdot ca- 
pacitance. The energy scale fku is, by comparison, only 
roughly known. From the fact that the barrier height Vq 
is approximately equal to E-p, we know that Xy ~ Af- 
For £, we can use the "device resolution" d, which is 
the distance between the surface metallic gates and the 
two-dimensional electron gas (2DEG) and is typically 
on the order of 100 nm. The fact that the approxi- 
mation £ ~ d should be accurate safely within a fac- 
tor of 2 can be surmised from calculations such as that 
of Davies and NixorO in which they show that the po- 
tential profile induced in a 2DEG by a narrow line gate 
has a half width at half maximum that is approximately 
equal to d.Ej'Eil In the AlGaAs/GaAs heterostructu*e£. 
of Waugh et ai, Crouch et ai, and Livermore et al.pu 
where d is fairly small, about 50 nm (approximately one 
Fermi wavelength), further circumstantial evidence for 
£ ~ d comes from the fact that the space between the 
gates that form the interdot barrier is about 100 nm (see 
Ref. 5). It follows that, for these experimental systems, 
hio is approximately 0.2_Ep- On the other hand, U is 
about 0.03-Ef, and, therefore, to within a factor of 2, 
2irU/huj ~ 1. For different systems in which the Fermi 
wavelength is stilLabout 50 nm but the gates are further 
from the 2DEG,E3 the ratio 2nU/fouj is presumably even 
larger. Consequently, we expect it to be quite generally 
true that the ratio U/W = 2ttII/Huj is greater than or 
approximately equal to 1. 

On the other hand, since, in the sorts of experimpru 
tal situations with which we are primarily concerned,cTu 
both W/Ep and U/Ep are much less than 1, we are jus- 
tified in linearizing the single-particle energy spectrum 
about the Fermi surface, taking E(j) = Ep + Hvp[k(j) — 
kp], where k(j) — ^/2mE(j)/h 2 . 

We must now calculate \B(j 2 , a 2 ;ji, ot\)\ when ji ^ j 2 . 
We avail ourselves of the exact, real solutions for the 
wavefunctions <ppj{x) in the presence of a parabolic po- 
tentiaMHl (P is the parity index, which we set equal to 
for symmetric wavefunctions and 1 for antisymmetric 
wavefunctions). Connecting these to the corresponding 
sinusoids, we find that, for x > \/2£ and Ldot ^> £, the 
eigenf unctions are approximately given by 



(~1) J 



dot 



ik(y Pj )(x-V2~o + Mypj)l (21) 



where ypj = y(Epj) and k(ypj) = y // 2mEp(j)/h 2 . The 
hard-wall boundary condition then demands that there 
be an integer n such that the quantity in brackets equals 
(2n + l)7r/2 when x = Ldot- 

As for the phase 7p(y) itself, it can be written in the 
following general form: 



ip{y) 



-l) p R(y) + D(y) 



(22) 



If the connection to the sinusoids is made using the lea cL 
ing large- a; forms for the parabolic cylinder functionsJij'cZI 
R(y) and D(y) are given by 



R( y ) = -arctan(e^) , 

D(y) = i [arg T(l/2 -iy)+y ln(47r V2£/A y )] + D , (23) 

where Do is independent of y. 

Returning to Eq. [l0| we find that, if we restrict the 
integral to x > \/2£, we have 

b' ~ (-i) q+i s j^M - D ^ cos ^(^) + r(vi)] 



2(k 2 - fci)idot 



B' ~ (-1) Q+1 C0S ^ D ^ Z D ( yi ^ sin [^(^ 2 ) + R (y^ (24) 

2(k 2 -ki)L do t 

where the bar of B 1 indicates that this is the term that 
moves electrons from dot a to dot a and we have antici- 
pated a continuum limit in replacing ypj and k(y P j) by 
Vj = yi E U)} and h = k{y 3 ). 

In the calculation of B 1 and B' , we have neglected the 
integration over the region \x\ < v2£- The results can 
therefore be expected to involve errors of order £/£dot 
when compared to the actual values of B and B. For B', 
this is not too much of a concern since, when both k 2 
and k\ approach the Fermi energy, the numerator of B' 
goes to ^fg and the denominator goes to zero. Thus, for 
non-infinitesimal g, if we restrict our wave- vectors to a 
range about the Fermi surface such that \ki — kp\ -C l/£ 
(in which case cos[D(y 2 ) — D(yi)] can be approximated 
by 1), corrections to B' should be relatively small. 

In contrast, the term B' is a bit more problematic, for 
its numerator goes to zero as {k 2 — fci)£. Consequently, 
near the Fermi surface this term is of the order of the er- 
ror, and to obtain a reliable result we must complete the 
integral numerically, using the parabolic cylinder func- 
tions in place of our sinusoids when |ar| < v2£- We 
then find that the form for B' approximates the mag- 
nitude of the actual value of B if, after approximating 
cos[R(y 2 ) + R(vi)] by 1, we replace [D(y 2 ) - D(y{)\ with 
K (y2 — yi), where k ~ 0.1 for g ~ 0.1 and n — > as g — > 0. 
Our conclusion is that 
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B ~ (-1) Q+1 sin ^ 2 Z 

2(/c 2 - fel)idot ' 

B ~ (-1) Q+1 sin [^(^) + 



2(fc 2 - fci)L do 



(25) 



We can now calculate the leading parts of the en- 
ergy shift A(p). The contribution from hopping between 
states on the same dot is given approximately by 



( u -f) 




{2wU\ 









Y, JO 



{V2 - yi) 3 



where the y^s are now measured relative to j/f (i-e., 
Vi — * 2/i + J/f)- The contribution from hopping between 
the dots obeys 



2 /-0 



dk x ( Al dk 2 



-A3 



sin 2 [i?(/i'i;Ffc2/^) + R(Hv-pki/ficj)] 
(fc 2 -fci) 2 
(1-P) 2 

fi,W F (fc2 - fcl) + C/(l - p) 

+ [p^-p}}, (27) 

where R(y) — R[u(Ef) +y] and the bracketed expression 
p — > — p stands for the quantity obtained by replacing p 
by — p in the previous term. In Eqs. ^6| an d |2^, ultravio- 
let cutoffs Y r and A r have been inserted in recognition of 
the fact that our formulas for the integrands break down 
at some distance Y r or A,, from the Fermi surface. 

(2) 

For the same-dot-hopping shift A c , the presence of 
such cutoffs is essentially irrelevant since we find this 
term to be effectively negligible no matter what the 
choice of Y r . In particular, even when the cutoffs are 
taken to infinity, this segment of the energy shift pro- 
duces a contribution to /(p) (recall Eq. |l^) that is 
bounded by the following formula: 



/c(p) 



< 



N ch p 2 f 2wU 



200 



(28) 



The real contribution is perhaps substantially smaller 
than the bound because the integrand of Eq. ^6|is system- 
atically too large for the infinitesimal-transmission states 
that correspond to \yi\ > 1. 

In any case, it is clear that the contribution to /(p) 
from same-dot hopping is essentially negligible. For 
A r ch < 2 and (2ttU '/ftw) <C 10, fc(p) is extremely small 
and essentially constant in g. Under such circumstances, 
it does not significantly affect even the quantitative re- 
sults. When {2 / KU/hu>) > 10, on the other hand, it can 
be relatively large. Nevertheless, it remains unimportant, 
for in this regime we can only obtain qualitatively good 



results for the value of A^(p), and, qualitatively, the up- 
ward shift of /(p) induced by same-dot hopping merely 
reinforces the effect from hopping between the dots. 

We now consider the cutoffs A r and their impact upon 
our understanding of the interdot-hopping result. Since 
the integrand in Eq. |2?j is reliably precise only when k\ 
and k 2 are within 1 /£ of &f , the ultraviolet cutoffs should 
be chosen such that A r ~ It follows that, to cap- 
ture with quantitative precision the leading behavior of 
the peak splitting for g <C 1, the set of wave- vectors 
within l/£ of the Fermi surface should encompass the 
range of energies in which the quantity R{E) is rapidly 
growing. Consequently, the set of wave-vectors must ex- 
tend at least to fcp + fcg, where E(k-p + ko) = Vq. From 
Eq. [l^, we see that the identity kg = (Vo~Ep)/hvF yields 



1 



2ttV2 



hi 



(29) 



If we require that fco£ < 1, we see that, for g <C 1, we 
must have \ hig\ < 2"K\[2. Thus, we have a lower bound 
on the values of g for which our approximations are reli- 
able. Fortunately, the lower bound is very small, and the 
requirement is only that g > 10~ 4 . 

( 2") 

We are now prepared to calculate A y T '(p). 
ter a switch to the dimensionless variables x r 
(-l) r hv F k r /U, Eq. [27] reduces to 



Af- 



A^p) 



N ch U 
4tt 2 
+ [p- 



f(xi,x 2 ) 

X-2 + X\ + 1 - p 



(30) 



where \r = fiv-pA r /U, the symbol = signifies equality 
modulo terms that are independent of p, and the quan- 
tity T{x\,x 2 ) is given by 

f{xi,x 2 ) = sm 2 [R(Ux 2 /huj) + R(-U xi/fuv)] . (31) 

To obtain a result with negligible dependence on the 
cutoffs Xn we must have Xr ^ 1< O n f ne other hand, 
to ensure that the answer is quantitatively reliable, we 
need A r < l/£ or, equivalently, Xr ^ hv-p/U^. Thus, as 
promised, we can only expect Eq. |3FJ to give quantita- 
tively reliable results for U <C Hvf/^; i.e., for 2ttU/Huj <C 
2tt\/2~, where 2nV2 ~ 9. 

Having dealt with the issue of the ultraviolet cutoffs, 
we can now go about the business of evaluating the right 
side of Eq. |o[ From the identity T(0, 0) = 5, it follows 
that the limit 2irU /frio — > yields the zero-width (£ = 0) 
linear-in-g equation for A^(p) that was pi|eyjausly de- 
rived via a transfer-Hamiltonian approach Jj&ffl In con- 
trast, in the limit 2-kU jhw — > 00, the energy shift given 
by Eq. 30 is independent of the interdot conductance for 
g a finite distance from both and 1. The constancy of 
the shift follows from the fact that, except when g equals 
or 1, T(xi,X2) is always 0.5 within the bounds of in- 
tegration, and the relevant parts of the energy shift are 
therefore the same as for (2-kU/Two) — and g = 0.5. 



G 



It should be re-emphasized, however, that such a result 
for the limit 2nU /hio — > oo can only be expected to be 
qualitatively correct. 

What happens when the barrier width £ is between 
and oo? By performing two partial integrations of the 
righthand side of Eq. [30| and dropping terms that go to 
zero as the cutoffs A r become infinite, we find that 



PH P ) = Nch9 



(l-p)ln(l-p) 



df(xi,0) 




h(p, xi,0) 



h(p,0,x 2 ) 



xi rX2 d 2 T(x ly x 2 ) 

dxi I dx 2 — s — ^ h{p,xi,x 2 ) 

^0 



dx\dx 2 



[P -> ~p] > 



(32) 



where h(p, xi,x 2 ) — {x 2 + x\ + 1 — p) h\(x 2 + x\ + 1 — p). 
The first term on the righthand side of Eq. |3^ is the zero- 
width result. The other terms, which go to zero in the 
limit £ — > 0, are the corrections due to a nonzero width. 

Numerical evaluations of Eq. [32] are plotted in Fig. 2(a) 
for several values of the parameter 2nU /ftto in addition 
to the analytically derived results for the limits of zero- 
width and infinite- width barriers. A curious feature of 
these curves is that the corrections to the zero-width 
behavior are antisymmetric about g = 0.5, a property 
that can be demonstrated by considering what happens 
to the integrands under the transformations <?<-> (1 — g) 
and X\*-^x 2 . Though the antisymmetry is suggestive, it 
must be remembered that / W (p) is only the leading term 
in a perturbative expansion about g = 0. The small pos- 
itive contribution to f(p) that comes from the formula 

(2) 

for the same-dot- hopping shift breaks this antisym- 
metry, and other higher-order corrections are likely to do 
the same. Nevertheless, some sort of rough antisymme- 
try about g = 0.5 is probably preserved, for, just as we 
find that, at small g, f(p) is enhanced by hopping con- 
nections to states with large transmission amplitudes, so 
we can expect that, at large g, f(p) is diminished by the 
fact that many of the occupied states from which one 
hops have transmission probabilities that are less than g. 

Such musings aside, we can gain further insight into 
the nature of our result for f^(p) by making a rough 
analytic approximation to the righthand side of Eq. |32j. 
To do this, it is best to return to Eqs. [l^ and |3(] and to 
derive the equivalent expression 



f (i) (p) 



1 , [ x * T( Xl ,x 2 ) 

ctoi / dx 2 - — — - 

Jo (x 2 +x 1 + l) 

I 

(x 2 + x-i + 1 — p)(x 2 + x x + 1 + p) ' 



(33) 
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Dimensionless interdot conductance g 
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Dimensionless interdot conductance g 

FIG. 2. (a) Plots of the leading g -> term of /, the 
fractional peak splitting, as a function of g, the dimension- 
less interdot channel conductance. Each curve corresponds 
to a different value of the quantity 2nU /hio (see legend on 
right). All curves are for the case of two interdot tunneling 
channels, iV c h = 2. The upward sloping solid line is the 
linear-in-g result that comes from considering an interdot 
barrier of effectively zero width (2nU /hu; = 0). The dashed 
and dot-dashed curves show the /-versus-g dependence for 
finite- width barriers with 2irU /hui taking on values from 
0.5 to 32. The horizontal solid line gives the leading term 
in the fractional peak splitting for an infinite- width barrier 
(2irU/huj — > oo). The curves can only be expected to be 
quantitatively accurate when 2-kU /huj <C 10. (b) /-versus-g 
results for the full domain of g when N c -^ = 2. The solid 
lines are the complete zero-width results in the weak and 
strong-coupling limits. These results contain both leading 
and subleading terms. 6 The plot for the leading zero-width 
term in the weak-coupling limit (g — > 0) is included as a 
dot-dashed curve. The small-dashed curve that extends 
from (<?,/) = (0,0) to (g,f) = (1,1) is an interpolating 
curve that is derived from the zero-width results. The 
long-dashed line is the 2ttU '/hu = 1 curve from Fig. 2(a). 
The stars, triangles, and squares symbolize different sets of 
experimental data, 2 ' 4 the squares being the most recent. 4 



We then postulate that, for small g, the magnitude of the 
f^\p) is largely determined by the portion of the inte- 
gral that corresponds to x 2 > xq, where xq = hv-pko/U 



7 



(recall fco from Eq. ^9|). For x 2 in this range, T(xi,X2) is 
on the order of 1 and therefore much larger than T(0, 0) 
when g -C 1. We label this high- energy portion of the 
double integral as fj^ p (fl)- Since, in this part of the inte- 
gral, T(xi,X2) varies relatively slowly between 0.15 and 
0.5, we approximate it by a constant Cf, where we take 
0.25. For xq > 1, we can drop the p's that appear 



Pi 

in the integrand of Eq. 



We then obtain 
N chP 2 



4ir 2 (x Q + 1) 



(34) 



From the identities / = /(l) and xq 
g)/g], we conclude that 



{hiu/2irU) ln[(l 



f (i) 

J hep 



iVch 

4tt 2 



2-KU/huj 



| In g | +2irU/hu; 



(35) 



This rough approximation to the leading behavior of 
the fractional peak splitting is only valid when xq > 1 
and fco < l/£. The condition on xo is necessary to 
justify dropping the p's from the integrand in the high- 
energy portion of the double integral. The condition on 
fco validates the approximate values for the magnitudes 
of \A(j2,OL2;ji,ai)\ that we use throughout our calcula- 
tion. The two restrictions together mean that the range 



of reliability for our approximation to is given by 



2irU 



< 



lng\< 27rV2. 



(36) 



For 2itU /flu ~ 1, our approximations are good for g be- 
tween a couple tenths and a few ten-thousandths. 

It is instructive to compare our result for f^J p with the 

zero-width fractional peak splitting, f^ - From Eqs. |l| 

and |35|, we see that, for 2irU/tkJ = 1, the ratio f^p/ fj;=o 
is about 0.6 when g — 0.1 and about 25 when g — 0.001. 
For very weak coupling (g -C 0.1), the correction to 
the £ = result is proportionately very large, and, as 
g — > 0, it dominates the behavior of the peak splitting. 
On the other hand, as g assumes more intermediate val- 
ues (g ~ 0.1), the results for £ = and £ ^ converge. 

A direct comparison of our results for with the full 
numerical results for which are plotted in Fig. 2(a), 
confirms that f^J p does indeed capture the essential /- 
versus-g behavior, particularly as 2irU/Hui becomes larger 
and the exponential enhancement of the tunneling am- 
plitudes becomes more important. The sharp increase in 
slope as g — > can now be understood as resulting from 
the fact that, in this limit, the high-energy portion of the 
peak splitting is proportional to (2nU/tuj)/\hig\. This 
proportionality also explains why the increase in slope as 
g — > becomes less dramatic as 2-nU jfrio decreases. The 
success of our rough analytic approximation supports the 
supposition that, for small g, a substantial portion of the 
peak splitting comes from tunneling into virtual states 
lying near or above the top of the barrier. 



Turning to Fig. 2(b), we now examine the significance 
of the calculated finite-width corrections in the context 
of what we know about the /-versus-g curve in the entire 
range from g = to g = 1, and we consider the implica^. 
tions of these corrections for the relevant experiments. Bel 
The long-dashed curve in Fig. 2(b) is the curve from 
Fig. 2(a) for the value 2nU/tiw = 1, which we believe 
to be appropriate for the cited experiments. The dot- 
dashed line is the leading-order- in-g, zero- width curve, 
which also appears in Fig. 2(a). The small-dashed curve 
in Fig. 2(b) is an interpolation for the entire zero-width 
/-versus-i? curve. This interpolation has been designed 
to match both the second-order-in-g calculation of the 
fractional peak splitting for weak coupling (g ~ 0) and 
also the two-term calculation for strong coupling {g ~ 1), 
which were obtained in Ref. 6 and are shown as solid 
curves in Fig. 2(b). The stars, triangles, and squares 
represent different sets of experimental data. 

For the particular value of 2irU jh\o that is illustrated 
in Fig. 2(b), we see that, although the finite-width cor- 
rection to / changes the answer by a large factor in the 
region of small g, the correction is small on an abso- 
lute scale. The difference between the dashed curve and 
the dot-dashed curve never exceeds 0.02 and therefore 
causes only a small correction to the overall shape of the 
/-versus-i? curve. Qualitatively, the correction due to the 
finite thickness of the barrier is quite similar to adding 
a small constant to / near g — and then decreasing 
the slope of the /-versus-g curve at small g. This qual- 
itative similarity follows from the fact that the region 
where / drops rapidly to zero, at very small g, is almost 
invisible in the plot. Consequently, the correction to the 
zero-width curve might be hard to distinguish from the 
effects of a small interdot capacitance, which have already 
been included in analyzing the data. We therefore con- 
clude that introduction of the finite thickness correction 
has little effect on the agreement between theory and the 
existing experimental data, for which 2ttU jh\o ~ 1. Nev- 
ertheless, such corrections may be important in future 
experiments. 



IV. CONCLUSION 



By developing a new approach to the coupled-dot prob- 
lem that relies upon the non-interacting, single-particle 
eigenstates of the full coupled-dot system, we solve for the 
leading correction to zero-width, wsak-jfiQupling results 
that were derived in previous work. mm The nonzero 
barrier width £ and finite barrier height Vq mean that the 
off-diagonal "hopping terms" vary exponentially with the 
energies of the states they connect. For a small interdot 
channel conductance (g <C 1), the resulting enhancement 
of tunneling to "high-energy" states above the barrier 
leads to an increase in the magnitude of the fractional 
peak splitting / observed at a given value of g. For a 
parabolic barrier, the magnitude of this increase grows 
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with the value of the ratio 2wU /tuo, where U is the inter- 
dot charging energy and u> is the frequency of the inverted 
parabolic well. Except in a very small region near g = 
where / behaves like (2irU/fkS)/\ \a.g\, the increase in / is 
accompanied by a decrease in the slope of the /-versus-i? 
curve. The effect upon the overall shape of the /-versus- 
g curve is not very substantial for experiments in which 
(2ttU /hu) ~ 1 but could be crucial in interpreting exper- 
iments involving significantly wider barriers. 

One might worry that the finite-width corrections to 
higher-order terms in the weak-coupling expansion could 
lead to a more dramatic alteration of the /-versus-<? 
curve. However, the corrections to such "large-g" terms 
should be muted by the fact that, as g increases, there 
is less difference between tunneling amplitudes between 
states at the Fermi energy and tunneling amplitudes be- 
tween a state at the Fermi energy and a state lying above 
the barrier. 

A more vital source of concern might be the treat- 
ment of the electron-electron interactions in the vicinity 
of the barrier. Clearly, the use of a sharp step function 
in the equation for n (recall Eq. ||) is an artifice. A more 
realistic model would account for the fact that, though 
electrons in and about the interdot channel still repel 
one another locally, their interactions with the rest of 
the electrons in the system are screened by the surface 
gates. 

Finally, one might wonder whether higher-order correc- 
tions to / preserve at least a rough antisymmetry about 
g ~ 0.5. We have seen that the leading small- g correc- 
tion, when directly extended to g = 1, changes sign and 
becomes negative for g > 0.5. Although a proper calcu- 
lation of the behavior at such large values of g requires 
consideration of higher-order diagrams, which we have 
neglected, we believe that the negativity of the correc- 
tion to / at large values of g is a generally right physical 
feature. When g is large and the reflection probability 
at the Fermi energy is therefore small, the energy depen- 
dence of the reflection amplitude, for £ ^ 0, should lead 
to a decrease in / as a result of the enhanced reflection 
coefficient for occupied states lying below the barrier. 
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